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This study presents the supersonic inductive Magnetohydrodynamic (SIMHD) 
generator and simulates its model through the use of COSMOL Multiphysics. 
This generator, as an inductive MHD generator, is suggested to address the 
problems associated with the conventional MHD generator. Since the proposed 
generator does not require a moving part to convert thermal energy to electrical 
energy, it is categorized as a direct energy convertor. The SIMHD generator 
consists of a converging-diverging duct and is divided into two sections at the 
diverging part by means of a diaphragm. Both of these sections are diverging, 
which makes it possible to obtain a high Mach number. In this regard, the 
performance of the SIMHD generator is studied by its mathematical modeling 
and numerical simulation using the finite element method. In addition, a 
sensitivity analysis is carried out on the design parameters. The results indicate 
that the proposed SIMHD duct design will increase the speed from 160 to 1300 
m/s, and decrease the temperature from 2100 to 1300 K. Moreover, 
considering a charge generation equation, the produced power for the resistive 
load 50 Q equals 25.3 kW. This generator also has the potential to be used in 
scramjets and ramjets. 


1. Introduction 


generates. In addition, since an electric discharge is 
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Since the beginning of the twentieth century, the 
applicability of Magnetohydrodynamic (MHD) Generators for 
the replacement of gas turbine-based energy conversion has 
been widely investigated [1-8]. Without the utilization of any 
mechanically moving part, the MHD power generation can 
convert thermal energy into electrical energy, and so this 
energy conversion system has attracted much attention [9- 
12]. These generators are not efficient unless the level of the 
charge concentration in the gas is increased to the 
conductivity level [7, 13, 14]. To achieve this end, the gas is 
highly heated and seeded with ionizing elements like alkaline 
metals [6, 9, 13, 15]. One of the main problems of these MHD 
generators is the high temperature that is required for 
ionizing the gas to operate [2, 9, 15]. Two other problems are 
also associated with these generators, which must be taken 
into consideration: (1) they require a high magnetic field 
(about 5 T) for which superconducting windings are needed 
[16], (2) the high-temperature plasma can destroy the 
electrodes which are in contact with it [16, 17]. To overcome 
these shortcomings, references [18-21] introduced a type of 
MHD generator. This model of MHD generator requires no 
external magnetic field as its conversion process proceeds on 
the basis of the induction which a non-static current 


responsible for the production of charge carriers, their 
concentration is not influenced by the temperature. 
Consequently, no electrode necessarily undergoes 
deterioration because of the magnetic coupling of the plasma 
and the electric load. The inductive gas-fed generator of this 
type of MHD generator, in addition to observing all the merits 
of the conventional MHD generator (e.g., the static conversion 
of energy and working with high temperatures), makes it 
possible to obviate the aforementioned shortcomings. In this 
model, an external electric field and a diaphragm that divides 
the field current into two parts split the charge carriers 
preceded by a set of electric discharges that ionizes the 
operating fluid. The unsteady nature of the electrical current, 
which results from the drag movement of charge carriers, 
leads to the induction of an electromotive force, a toroidal coil 
that is put around the pipes. To increase the gas speed and 
decrease the distance between the electrodes, which results 
in a lower potential difference, a converging-diverging duct 
was utilized in this study. To identify whether a desirable 
level of energy can be transmitted to the electric load, some 
studies which involve this process have been carried out [18- 
21]. In this paper, the Finite Element Analysis (FEM) has been 
utilized for the investigation of Supersonic Inductive MHD 
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generators. In addition, to make the simultaneous 
consideration of both dynamic and_ electromagnetic 
dimensions, a multi-physics approach has been employed. In 
this regard, the physical phenomenon of the problem is 
described in Section 2. The numerical model is presented in 
Section 3. The model results are discussed in Section 4. The 
conclusion of the study is presented in Section 5. 


2. Physical phenomena description 

Like an inductive MHD generator, the supersonic 
inductive MHD (SIMHD) generator can be categorized as 
direct energy conversion since it directly converts the energy 
from fossil fuel into electricity. Like all the other direct energy 
conversion processes, this model of generators is capable of 
converting thermal energy into electricity without using any 
moving parts. In comparison with the conventional MHD 
generator, there is no need for an external magnetic field in 
this model. While due to the principle of inductivity, the 
conversion of energy is still possible. This generator operates 
based on simple principles. Initially, a high-speed gas enters 
the duct. Here, through a pulsating electric discharge, the gas 
is ionized. An external electric field is utilized to split the 
charges of the different signs. A force F = QE, which is 
parallel to the direction of E and separates the charges with 
different signs in the direction of the movement of the fluid, is 
produced on each charge Q of the fluid by the movement of 
the ionized fluid in an electric field Ë. The flow is divided into 
two charge currents by a diaphragm, one with an excess of 
positive charge and the other with an excess of negative 
charge. An electromotive force is induced in toroidal coils put 
around a highly permeable magnetic core which is placed in 
the cavity and around the duct by the generated current. 
Figure 1 displays a representational scheme of the SIMHD 
generator, which consists of three main sections. The first 
section is responsible for the ionization of the inlet fluid. In 
the second section, the charges are separated. Finally, in the 
third section, the energy is transmitted to the electric load. 
There are different methods that can be utilized for the 
ionization of the fluid. In this study, it is assumed that the 
electrodes, which are powered by a high-voltage generator to 
produce a pulsating electric discharge, ionize the fluid. 
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Figure 1. The SIMHD generator functional scheme 
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To reach the appropriate electric current, there is a need for 
the generation of an adequate number of charge carriers 
during the discharge process. The electrode must be located 
in the initial part of the duct, connected to the pulse generator, 
which produces a high voltage. It also needs to be in contact 
with the ionized fluid. It must also be noted that the pulse 
generator needs to be appropriately set in order to select the 
discharge generation optimal parameters. A DC high-voltage 
generator externally powers the two plates of a capacitor 
which are placed in the second section of the generator to 
separate the charges. To avoid the blockage of the charge 
carriers inside the device, the external electric field should be 
well-calibrated and fine-tuned. That is, the capacitor plates 
must be powered by a voltage that is sufficiently high to 
separate the positive and negative charges and not so high to 
block the charges in the internal side of the apparatus. Finally, 
two currents of fluid, one with an excess of positive charge 
and one with an excess of negative charge, enter the third 
section of the apparatus. The diaphragm is placed in the third 
section of the generator, and as soon as the fluid currents 
enter the section, they are divided into two asymmetric flows, 
which are electrically unbalanced. Therefore, a ring of electric 
current is created around the diaphragm (Figure 1). In this 
section, some of the energy of the fluid is carried to the 
external load, and it is here that the energy is converted. The 
ring of the generated electric current varies with time, and as 
a result, an electromotive force is induced in the toroidal coil 
around a core that has a high magnetic permeability. This core 
is around the generator in the cavity. The charges are then 
neutralized as the two currents combine downstream of the 
third section. It can be observed that an external electric field 
separates the negative and positive charges, and there is no 
need for the generation of any external magnetic field. What 
makes energy conversion possible is the existence of the 
inductive process. As an electric discharge produces the 
electric charge, the gas can be ionized at low temperatures if 
the pulsing voltage is sufficiently high, and since there is the 
magnetic coupling of the plasma and the electric load, no 
electrode will be destroyed. 
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3. Model description 

The system of coupled (nonlinear) equations is used to 
describe the mathematical time-dependent model of plasma 
flow in an electric field. This system is composed of Poisson's 
equation for electric potential, the Navier-Stokes equation 
and the conservation energy equation for fluid, and the charge 
transport equation for electric charge density. 


3.1 Electrostatic 
The Poisson's equation which results in the electric 
potential V is as follows: 
vVy=-+ (1) 
E0 
where, q (sı unit: Se is the space charge density , and £o 
is the dielectric permittivity of the free space. 


The electric field intensity E defines the electric potential, as 
in Equation (2): 


E = -W (2) 


3.2 Fluid dynamic 

The Navier-Stokes equations depict that part of the 
problem which is related to the fluid dynamic of the gas flow 
with the steady state turbulence. Turbulence modeling has a 
crucial role to play in the calculations of hypersonic flows. A 
turbulence model must always be included in the flow 
modeling since turbulent flows exist in many hypersonic 
applications. According to the one differential equation, 
Balswin and Barth have presented a new class of turbulence 
models. The most recent development of these models is 
represented by the Spalart-Allmaras (S — A) model. The 
classical Prandtl approach in which the turbulent kinetic 
energy equation is used to arrive at the model equation ,is not 
adhered to in this new class. To create the equation in the S — 
A model, Galilean invariance, selective dependence on the 
molecular viscosity, dimensional analysis empiricism and 
argumentation are employed. But the equation in Baldwin- 
Barth model is obtained from k— e model [22, 23]. The 
velocity, pressure, and temperature fields are obtained by the 
S — A model all the way down to the wall. That is why this 
model is referred to as a Low-Reynolds number model. 
Different applications of the S — A model has proven to be 
successful [24]. To be solved, this model only requires one 
transport equation of the quantity 7 (or vt), which equals the 
eddy viscosity v, far away from the walls. The reason behind 
the construction of this transport equation has been the 
empirical reproduction of the flows with great intricacy [25]. 


pCi: V)i = V: |-p? + (ut wr)(Vi + (Vi) - (u + 
ur)(V- a)i] +Ë (3) 


where, F (sı unit: Nia] is the body force and is obtained by 
the following equation: 


F =] xB- qvV (4) 


The Navier-Stokes equation does not take the body forces into 
equation. With the left circuit open, the magnetic field 
reaction is near zero (Jx B= 0) in the case investigated. In 
addition, the body force (—qVV) can be ignored since the 
concentration of charge q in the fluid is very low. The 
relationship between the shear rate and the shear stresses in 
a fluid is described by the dynamic viscosity u(SI unit: Pa: s). 
The Sutherland's law is used to compute the coefficients of 
viscosity [25, 26]. The relationship between the ideal fluid's 
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dynamic viscosity and its total temperature is described by 
this law based on: 


3/ 
T 2 TureftSu 
T+Sy 


u= Hrer ( (5) 


Turef 


where Href is the dynamic viscosity at reference temperature, 
Tyrer (SI unit:K) is the reference temperature, and 
S„ (SI unit: K) is the Sutherland constant. Table 1 presents 
Sutherland’s law parameters for dynamic viscosity for 
different gases. In this study, the properties of the air were 
utilized to solve the equations of the fluid dynamic. 


Table 1. Sutherland’s law parameters for dynamic viscosity 
(adopted from COMSOL software's documents) 


Hre, Ture 5 
cS (Paes) CK). (K) 
Air 1.716 x 1075 273 111 
co, 1.370 x 1075 273 222 
co 1.657 x 107 273 136 
N3 1.663 x 1075 273 107 
0; 1.919 x 107-5 273 139 
Steam 1.12 x 1075 350 1064 


A partial differential equation is employed to calculate the 
eddy viscosity in the S — A model. The parameters used in the 
S-A model is presented in Table 2. An intermediate variable 
vt is used to calculate the viscosity ur through the following 
equation. 


br = pV fo X) (6) 
where y is the ratio determined by Equation (7): 


vt 


y (7) 


and f,,1 is a damping function calculated by Equation (8): 


foi = 


3 


(8) 


X 
Ot, 


ae ce rors (9) 


n2 
T.V = Cy St" — Cui fy (=) # av ‘(vy + vt) Vv) + 


ca Yyt . Yyt yt = nutilde (7) e 
Cw = Zz + aa = 


Table 2. S — A model’s parameters (adopted from COMSOL 
software's documents) 


Constant Value 
Cpi 0.1355 
Cho 0.622 
Ĉi 71 
Oz 2 J3 
Caz 0.3 
Cas 2 
Ky 0.41 
Crot 2 
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1 
6 


= 1+C3 
haal) (12) 
g =r — Cm(rf -r) (13) 
; vt 
r= mine’ 10) (14) 


t 
St = max (Q + Croemin(0, S — 2) + aga 032) (5) 


2=V20:0 (16) 


S=V28:8 (17) 
0 = (vù - vÙ") (18) 
Š=: (yù + (Vi)") (19) 


The mass conservation is denoted by Equation (20): 
V: (pu) =0 (20) 


The conservation of energy can also be described by Equation 
(21): 


PCy VT =V: (kVT) + Qp (21) 


where cp (SI unit fog . K) is the heat capacitance of the gas 


at constant pressure, and Q,(SI unit: W/ 5) includes the 
heat sources, determined by Equation (22). 


Q=]: B (22) 

If the secondary circuit is open, there will be no transfer of 

energy. The heat sources can be neglected in the evaluation of 

electrohydrodynamic flow under these circumstances. 

For thermal conductivity the Sutherland’s law can also be 

formulated as follows [26]: 
ire 

T 2 Tere tSk 
ae (23) 


k = krer ( 


where kref (SI unit: Wr ; K) is the thermal conductivity at 
reference temperature, Ty;er (SI unit: K) is the reference 


temperature, and Sx (SI unit: K) is the Sutherland constant 
which is mentioned in Table 3. 


T ref 


Table 3. Sutherland’s law parameters for thermal 
conductivity (adopted from COMSOL software's document) 


kref 


Tkref Sk 
ca: —) (K) (K) 
m:K 
Air 0.0241 273 194 
co, 0.0146 273 1800 
co 0.0232 273 180 
N3 0.0242 273 150 
0, 0.0244 273 120 
Steam 0.0181 300 2200 


Finally, the equation of state is as follows: 
p = pRT (24) 


3.3 Charge transport 

Usually, the charge carriers which are present in the fluid 
are electrons and ions. The complex behavior of electrons in 
liquids or gases may be due to the collisions of ions with 
neutral molecules and other ions, photoionization, etc. the 
calculation of the velocity of the ions is a daunting task, and it 
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has a great dependence on the local electric field [27, 28]. If 
the electronegative impurities or the molecules present in the 
fuel capture the electrons, then the only charge carriers in the 
fluid are the positive and negative ions. Under such 
circumstances, it is possible to regard these ions as spheres 
which are moving in a continuum and estimate their behavior 
based on the laws of an ideal gas of particles in the volume 
which is filled with the fluid [28]. Three factors affect the 
electric current in the drifting zone: conduction (motion of 
charges under electric field relative to entire flow), 
convection (transport of charges with flow), and diffusion [28, 
29]. As a result, the current density J is calculated by Equation 
(25): 


J = ugËq + qu—DVq (25) 


Combining equation (2) with equation (25) leads to the 
following equation: 


J = q — DYq — urq VÝ (26) 


where, Hg (sı unit: m? hy . 5) is the mobility of the charge in 
an electric field, ŭ is the flow's velocity vector, and 
D (sı unit: miz) is the diffusion coefficient of charges. By 


making a comparison between the peripheral particles' 
diffusion velocity and the velocity of these particles which 
results from their interaction with the distribution of the 
charge, the value of the diffusion coefficient is computed and 
by increasing this value, the effect of the charges repulsion has 
been accounted for. Taking into consideration the charge 
carriers' mobility, the diffusion would be D = 2.1208 x 


1075 m? |. instead of D=5x10-5 m? |< which is the 
diffusion without charges interaction [19, 21, 29]. The 
equation for current density based on the current continuity 
condition is as follows: 


aq = q For charge generation region 
—+V-'J= (27) 
0 For other regions 


where, ġ (sı unit: Chn . J is the charge generation rate. 


a4 V: (qu — DVq — ueq V7) = 

q For charge generation region 

(28) 

0 For other regions 

3.4 Charge generation region 
A gasifier or a combustion chamber feeds the duct with a 

gas [4, 30]. The literature indicates that the concentration of 
the charge carriers which results from the discharge is around 


1014 "ey n3 [31, 32]. Assuming that in the process of 


discharge 10 cm? is engaged, the overall number of charges is 
approximately 2 x 1015. n,-e = 1.6 x 1074C and it will 
then be the charge of the same sign. A time varying flow of 
positive and negative charge carriers is injected to model the 
generation of the charges by the electrode. A space charge 
density q is used in the electric discharge region by applying 
the reaction rate parameter to the equation that follows: 


__ 2 (t-t)? | (t-t)? q = (ttn)? 
2162 pC 262 t 262 T t 262 ) (29) 


where, T 3.) is the charge generation 


function, 7 and 6 are the two parameters for adjusting this 
equation with charge density distribution, and ty, t2,...,t, are 
the times when there is the maximum concentration of the 
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charges. In this context, Figure 2 displays the charge 
generation function and the time variation of charge density. 


x104 
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T T 


N 
a 


a 
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Charge Generation Function(C/m?.s) 
= N 


o 
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t(s) x10 


x10183 


e 


Charge density(n /cm*) 


(b) 


Figure 2. Time variation of (a) charge density function and 
(b) the charge generation function (7 = 30,6 = 2.8 x 1074) 


3.5 Energy conversion process 

The functioning of the third section of the duct, where the 
physical process of energy conversion takes place and the 
electrical energy is transferred from the fluid electric current 
(primary winding) to the secondary winding by means of the 
inductive linkage, is the same as the functioning of a 
conventional transformer [33]. To calculate the force of the 
inductive electromotive force, an equivalent circuit in which 
the primary winding is equal to the current ring which flows 
along two parallel branches is assumed. It is also assumed 
that the secondary winding is comprised of a number of 
windings around two ferromagnetic rings which are put 
around the aforementioned branches (Figure 3). A varying 
electromotive force is induced in the secondary winding by a 
varying magnetic flux which is created by a varying current in 
the primary winding. An electric current flow in the 
secondary winding provided that a load is linked to it. In 
addition, an electric current will be transported to the electric 
load from the primary circuit through the transformer. An 
armature reaction is created by this electric current. This 
reaction reduces the speed of the charge carriers in the fluid 
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and leads to the expansion of the gas. The magnetic flux 
(SI unit: Wb) is obtained by the following equation: 


(t) 
t= (30) 


where Rm(SI unit: H/m) is the magnetic reluctance. The 
reluctance of a magnetically uniform magnetic circuit element 
can be computed as: 

L 
Rm = a (31) 
where (SI unit:m) is the length of the element, 
Lm (sı unit: N/ 42) is the permeability of the material and 


A(SI unit:m?) is the cross-section area of the circuit. The 
permeability of the material is determined using Equation 
(32): 


Um = Emo * Umr (32) 


where [mo = 4m X 1077 Ne is the permeability of free 


space, and Hm, is the ratio of the permeability of a specific 
medium to the permeability of free space that is assumed 
equal to Um, = 4 x 104. 


2RmitRm2 


Rm = Rmo + (33) 
2I 

$: = 2Rm,1+Rm2+2Rmc (34) 
doc 

emf = Ns (35) 


where N, is the number of coils around the magnetic core. 


Rai (36) 


Power = vn, — == 
rms Toa 


where V,ms is the root mean square voltage and Rg, and Rg, 
are the external load resistive and the coil resistive 
respectively. Ifassumed that Rgı = Rg ¿ equation (36) will be: 


vms 
Power = Tre (37) 


Yims = |+ p" emf)?dt (38) 


where tm = ti+1 — t; is the time between two consecutive 
electric discharges. 


Rm.2 
Lc, L2 


Figure 3. The equivalent magnetic circuit and 3D magnetic 
circuit scheme 
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4. Results and discussion 
4.1 Modeling results 

A 2D configuration model has been utilized in the present 

study. A converging-diverging shape is applied to the volume. 
The duct has a height of 110 mm, a length of 234mm and a 
throat height of 40 mm. Figure 4 displays the boundaries, 
mesh, numerical simulation space, and subdomains of the 
model. The COMSOL Multiphysics 5.0 was employed to run 
the model by using the Finite Element Analysis (FEM). Three 
separate modules were used as follows: 

e Electrostatics to solve the equations related to the 
electrostatic section; 

e Fluid Flow (High Mach Number Flow (Turbulence Flow, 
Spalart - Allmaras)) to solve the Navier-Stoke's and energy 
equations; 

e Transport of Diluted Species module that addresses the 
convection, diffusion, and migration in electric field 
aspects. 

To solve the charge transport equation (Equation 28), the 

Diluted Species module was utilized. Since this module was 

used and the dependent variable of this module is the 

concentration of the components (c;), where its unit 


is (ee a) the concentration of the component is 


considered to be the same as the charge density of the 
components (positive ions and electrons) q;. Following this 


assumption, ae can also be considered the same as 


(c m3) 

A stationary study was carried out to solve the first two 
modules. The results of this study were utilized in a transient 
study to analyze the third module. To solve the equations of 
the charge transport in which transient interactive solvers are 
employed, the velocity and the electrical fields are computed 
and utilized as source terms. In the stationary study, the mesh 
in Figure 4 was used, which is appropriate for solving fluid 
equations. But in the transient study, the mesh in Figure 5 was 
used. The mesh depicted in Figure 4 is comprised of nearly 
130000 elements which are more concentrated near the 
walls. 


Figure 4. the boundaries, mesh, numerical simulation space, 
and subdomains of the model. (x and y axises are in mm) 


An injection of charge carriers was done to simulate the 
electric discharge, which is shown by the circle. In this region, 
Equation 29 was inserted as the entry. The semi-circle around 
the discharge region was used as a virtual boundary to refine 
meshing. The following boundary conditions were applied to 
the numerical model. A no-slip function was employed in the 
inner walls of the duct in the fluid dynamic study. In the inlet 
section, the subsonic condition was applied to the inlet witha 
velocity of M = 0.2 and a total pressure of 5 atm. In the outlet 
section, the hybrid condition was applied to the fluid with a 
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total pressure of 1 atm. A positive 30 kV DC voltage, which 
was constant, was applied to the external plates for 
electrostatic analysis. The mesh in Figure 5 consists of 
110000 elements which are more concentrated where the 
charges move. This mesh was applied to the Transport of 
Diluted Species model where a zero diffusive flux condition is 
applied to all the boundaries but the outlet section in which 
convective flux diffusion is applied. Table 4 shows the 
modeling parameter values employed in FEM modeling. 


° 20 40 60 80 100 120 140 160 180 200 220 


Figure 5. The mesh used in the charge transport equation 


Table 4. The modeling parameter values used in FEM 
modeling 


Parameter Value 


Relative dielectric permittivity of air, ¢,. 1 


2 
Charge diffusion coefficient, D 2.12 x 1074 M/s 


pant 
Ion mobility coefficient, 4g 1.8 x 10°-4™ ly -s 


k 
Density of the fluid at the inlet, Pfiuia 1.23 Ilai 
Number of the coils around the magnetic 8000 
core, Ns 
Time between two consecutive electric 0.0023 s 
discharges, tm 
Height of the cavity, H 8mm 
Mach number of fluid at the inlet, Min 0.2 
Inlet total temperature 2000 K 
Inlet total pressure 5 atm 
Outlet total pressure latm 
The difference in voltage between the 30kV 


capacitor plates 


As it is shown in Figure 6, the highest potential of the electric 
potential map is at the surface of the top capacitor armature, 
while at the bottom capacitor it has a decreasing magnitude. 
The magnitude of the gas velocity inside the channel is 
displayed in the velocity field map (Figure 7). As it can be seen 
in Figure 7, the oblique shocks can be observed in the regions 
where the velocity is higher than 1 Mach. These regions are 
shown specifically in Figure 8. The temperature and pressure 
fields are shown in Figure 9 and Figure 10, respectively. The 
oblique shocks can also be seen in these two figures. The 
movement of the charge carriers inside of the duct and the 
magnitude of the electric current density is shown at the 
5 x 10-*th, 1 x 10~3th, and 2 x 10-3th seconds in Figure 11, 
Figure 12, and Figure 13, respectively. In Figure 14, the 
electric field lines, streamlines, and the direction of the charge 
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Pressure (atm) 


carriers' movement can be simultaneously observed. In this | á 6.2377 


figure, the effect of fluid's movement and electric force on the 100 ‘Je 
movement of charge carriers is also shown. 90) | 
Electric potential (V) a | 
A 3x10! ai | 
30000 e 
50| l 
25000 a | 
| 30) | 
20000 20 | 
10 1 
15000 o 1 1 


° 20 40 60 80 100 120 140 160 180 200 220 V 0.5872 


u 


= 
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Figure 10. The pressure field in the (atm) 
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110 4 133.37 
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Figure 6. The distribution of electrical potential in the duct ji haó 
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110 A 1358.2 50 | s 
100 | 40 | 
90 | 91200 = | 40 
20 
a | 1000 10 | 20 
70 | o | 
60 | [1800 o 20 40 60 80 100 120 140 160 180 200 220 V 6.9062x10% 
50 4 
600 
40 C 
30 | faco Figure 11. The electric current density ( [m : 5) at0.0005s 
20 { 
10 po 
o 1 
ò 20 40 60 80 100 120 140 160 180 200 220 Y 7.7444x10" Time=0.001 s sqrt(tds.tflux_qlx*2+tds.tflux_qly*2)+sqrt(tds.tflux_q2x*2+tds.tflux_q2y*2) (mol/(m**s)) od 
110 A 706.94 
100 700 
. . a . 90 
600 
Figure 7. The velocity field in the duct (m/s) 7 | 
70 500 
. Contour: Mach number 60 400 
110} A2 50 
n2 40 300 
100: 
90} ( 7 200 
1.6667 20 
80} 10 100 
70} H1.3333 6 
60! ~ | o 20 40 60 80 100 120 140 160 180 200 220 Y 1.0255x10™® 
50} Ha 
sl x F F C 
20! 0.6667 Figure 12. The electric current density ( [ua i 5) at0.001 s 
20) 
0.3333 
10; 
oj +o 
o 20 40 60 80 100 120 140 160 180 200 220 Y 0.3333 Time=0.002 s sqrt(tds.tflux_qlx*2+tds.tflux_qly*2)+sqrt(tds.tflux_q2x*2+tds.tflux_q2y~2) (mol/(m?*s)) a 
110) i i ` i - - | & 10.836 
Figure 8. The mach number contour in the duct set | fio 
" | 9 
80 
8 
70 z 
Temperature (K) ° 
` j ` s 7 7 p à A 2126.1 “ 6 
50 H 
|i 2100 5 
40 { 
t 2000 4 
30 | 
f 1900 3 
} 1800 a i E3 
j 1700 ao 1 
i 1600 z 
[o] 20 40 60 80 100 120 140 160 180 200 220 V 6.9617x10" 
if 1500 
1400 
ene Figure 13. The electric current density Ci 2 at 0.002 s 
1200 ms 
1100 
V 1087 


Figure 9. The fluid temperature field in the duct (K) 
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Time=0.001 s sqrt(tds.tflux_q1x*2+tds.tflux_qly*2)+sqrt(tds.tflux_q2x~2+tds.tflux_q2y~2) (mol/(m**s)) Velocity field Electric field 
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Figure 14. Electric field lines (red), streamlines (black), and 
the direction of the charge carriers' movement (the colored 
region, the unit of which is (C/(m2:s))) 


An equivalent circuit, in which the motion of the electric 
charge in the fluid creates the primary winding and the 8000 
coils, which are distributed over the branches of the magnetic 
circuit, constitute the secondary winding, is assumed so as to 
compute the inductive electromotive force. As shown in 
Figure 15, a voltage on the tops of the coil is arrived at by 
assuming the secondary winding is an open circuit. The 
assumed maximal power output has been computed by 
considering an adjusted load. 

A theoretical maximum power near 25.3 kW is achieved by 
assuming Eq (36). By considering Rz; = Rg c as mentioned 
before, the maximum power has been computed with the 
assumption of an adapted load. Under these assumptions, 
several factors reduce the functioning of the device, for 
example, if the external resistive load or the time interval 
between two successive electric discharges increase. The 
amount of power generated by changing the resistive load 
and the time interval between two successive electric 
discharges is shown in Figure 16. Here tm represents the 
minimum time interval of interference between the two 
electric discharges. The decrease in the performance of the 
device is due to the increase of the external resistive load 
from 51 to 800 and the extension of the time interval 
between two electric discharges from tm to 2.00 tm is shown 
in Figure 16. 
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Figure 15. The inductive EMF on the heads of the open circuit 


4.2 Sensitivity analysis 

Taking some design parameters, such as the inlet 
temperature ofthe gas and the geometry of the generator into 
consideration, a sensitivity analysis was carried out. The 
functioning of the voltage of the coil's head and the theoretical 
maximum power load is calculated to see the effect of these 
parameters. The profile width of the diaphragm and the 
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internal cavity of the diaphragm are the major dimension 
parameters that are taken into consideration. The working 
inlet temperatures have been changed at four points from 
1800 to 2500 K to see how they affect the voltage of the head 
of the coil and the theoretical maximum power load. Figure 
17 displays the voltage of the open circuit, which is obtained 
from changing the temperature of the working fluid and the 
height of the internal cavity of the diaphragm. As the 
temperature rises, the velocity of the charges increases, and 
as a result, the inducted gas voltage rises. However, it takes a 
shorter time period for the charge clouds to disperse; 
therefore, the charge clouds stay more concentrated. That 
section of the core which is inside the cavity becomes a 
bottleneck for the magnetic circuit and restrains the 
magnitude of the magnetic flux. 
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Figure 16. The trend of theoretical maximum generation 
power (kW) at different resistive loads (Q) and time intervals 
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Figure 17. The voltage (V) of the open circuit obtained from 
varying the temperature (K) of the working fluid and the 
height of the internal cavity of the diaphragm 


5. Conclusion 

In the inductive MHD generator, there is no need for an 
external magnetic field, and it can perform at lower working 
temperatures in comparison with conventional MHD 
generators. The innovation of the SIMHD, in comparison with 
the inductive MHD generator, is that it uses a converging- 
diverging duct which increases the velocity, which in turn 
increases the voltage and generation power. This increased 
velocity will also result in a decrease in temperature in the 
second section of the duct, where the charge carriers are 
separated and cannot combine together. In the duct designed 
in this study, since the distance between the capacitor plates 
is less than that of the inductive MHD generator, a lower 
potential difference is required to separate the charge 
carriers. The Finite Element Analysis (FEM) and a multi- 
physics method which makes the simultaneous investigation 
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of the aspects of high-Mach number flow, electrostatic and 
charge transport possible, were utilized to investigate the 
SIMHD generator in this study. COMSOL Multiphysics was 
used to model the present study. The results showed that the 
proposed SIMHD duct design increases the speed from 160 to 
1300 m/s, while decreasing the temperature from 2100 to 
1300 K. This can clearly have a positive effect on the 
performance and lifetime of the proposed generator. 
Moreover, the produced power for the resistive load 50 Q 
equals 25.3 kW, which is a desirable and acceptable number. 
The present study also opens the venue for further research 
in this area. Future studies can work on the applicability of 
this generator in scramjets and ramjets since they make it 
possible to reach a speed of more than 1 Mach. 
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